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Abstract 

TARCER is an implementation of the recurrence algorithm of O.V. Tarasov for the 
reduction of two-loop propagator integrals with arbitrary masses to a small set of 
basis integrals. The tensor integral reduction scheme is adapted to moment integrals 
Ctf ' emerging in operator matrix element calculations. 
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1 Program Summary 

Title of program: TARCER 
Version number: 1.0 



Available at [http : // www . mert ig . com/tar cer| or 



http : //www . physik . uni-leipzig . de/TET/ tarcer 



Programming Language: Mathematica 3.0 
Platform: Any platform supporting Mathematica 3.0 
Keywords: Feynman integrals, two-loop corrections 

Nature of physical problem: Reduction of Feynman integrals in pertubative 
two-loop calculations 

Method of solution: The system of recurrence relations given by Tarasov [1] 
Limitations: Rank of integrals 



2 Introduction 



TARCER reduces two-loop propagator integrals with arbitrary masses to sim- 
pler basis integrals using the reduction algorithm proposed by Tarasov [1,2]. 
The reduction of scalar products in the numerator of the integral is extended 
to provide for the presence of an additional external vector that is lightlike. 
For the reduction of scalar integrals TARCER contains the complete set of 
recurrence relations given in [1] and some additions for particular parameter 
configurations. In some cases the number of basic integrals is reduced. Further 
additions may easily be added. 

Due to a vast number of interrelations between the integrals considered it is not 
immediately obvious how to extract a set of recurrence relations that reduce 
the complexity of the integrals at each step such that one finally arrives at only 
a small set of basic integrals^. This was achieved by Tarasov. The resulting 

3 The fact that in principle all integrals considered must be expressible in terms of 
a finite set of basic integrals may be seen as follows [3]: Interpret the corresponding 
Feynman parameter integral representations as integrals over projective differential 
forms. All forms in question exist on the same differential manifold. De Rham co- 
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relations are in some cases quite involved. Therefore transcribing them into a 
program is an error prone process. But, whereas the extraction of adequate 
relations proves to be quite difficult they can be checked in most cases much 
more easily. For four or less propagators this may be done via a Mellin-Barnes 
representation of the integral. In the present program this is automatized in 
order to minimize input errors. 

The emphasis of TARCER is not so much on the speed of evaluation. Our 
main purpose was to use TARCER to check complicated moment integrals 
emerging from operator matrix element calculations for small individual mo- 
ments (to < 5). An example of this type of integral is given in section 7. 
Conversely, those integrals, calculated independently by other means, served 
as a check for TARCER. 

The massless case was checked using a FORM program [4]. The reader may 
also want to compare this program with other existing fully analytic two-loop 
programs like [5], [6] and [7,8] or the original implementation in [1]. 



3 Notation 

TARCER reduces the following general type of integrals to basic integrals. 
The Mathematica notation, making all dependencies explicit, is listed first: 

TFl[d, p 2 , Ap, {a, b}, {u, v, r, s, t}, {{v u mj, {{u 2 , m 2 }, . . . , {z/ 5 , m 5 }}] = 
1 r r d d k x d d k 2 (Ah) a {Ak 2 ) b (kp u {k 2 2 ) v {phY {pk 2 ) s (hktf 
n d J J [kf - to?]" 1 [kl - ml]"* [kl - m§]"s [kj - ml] u * [kl - ml]^ ' ^ ' 

with the abbreviations k^ = k±—p, k± = k 2 —p and k§ = k\ — k 2 . The exponents 
a, . . . ,t and the indices vi, . . . , v 5 are assumed to be nonnegative integers. A 
denotes a lightlike vector with A 2 = 0. 

If some of the subsets {a, b} or {u, v, r, s, t} of exponents vanish we have the 
following reduced notations: 

TFl[d, p 2 , Ap, {0, 0}, {u, v, r, s, t}, {{v u mj, {u 2 , m 2 }, . . . , {u B , m 5 }}] = 
TFl[d,p 2 ,{u,v,r, s,t},{{z/ 1 ,m 1 },{z/ 2 ,m 2 },...,{z/5,m 5 }}] , (2) 

TFl[d, p 2 , Ap, {a, b}, {0, 0, 0, 0, 0}, {{v u mj, {u 2 , m 2 }, . . . , {u 5 , m 5 }}] = 
TFl[d,p 2 , Ap, {a,b}, {{z/^mi}, {z/ 2 ,m 2 }, . . . , {z/ 5 ,m 5 }}] , (3) 

homology then implies that there exists only a finite number of inequivalent forms 
on this manifold. This in turn implies the above statement. 
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TFl[d, p 2 , Ap, {0, 0}, {0, 0, 0, 0, 0}, {{u u mj, {u 2 , m 2 }, . . . , {u s , m 5 }}] = 
TFl[d,p 2 ,{{z/ 1 ,m 1 },{z/ 2 ,m 2 },...,{z/ 5 ,ni 5 }}] . (4) 



Furthermore, if a mass vanishes the argument {uj, 0} of TFI may be replaced 
by the index Vj alone. For example the fully massless integral with five proa- 
gators may be entered as TFl[d, p 2 , {1, 1, 1, 1, 1}]. 

In the course of evaluation all scalar products in the numerator will be elim- 
inated. If no more scalar products are present we keep the notation in line 
with [1]: 



TFl[d, p 2 , {{v u mj, {u 2 , m 2 }, . . . , {u s , m 5 }}] = 
(d) 1 f f d d kid d k 2 

K,u m u 5 = — d j] y k 2 _ m 2^ Vl [£2 _ m 2y 2 ... [ fc 2 _ m 2], 5 > < '< 



TVl[d,p 2 , {{z/^mi}, {z/ 2 ,m 2 }, {z/ 3 ,m 3 }, {z/ 4 ,m 4 }}] = 
V ( d ) 1 /■ / d d hd d k 2 



(6) 



= _ / / a «i fl «2 m 

^1^3 / / TL2 ™21i/i TL2 ™21ia> f7„2 ™21;y, ' V/ 



TJl[d,p 2 , {{i/ijitti}, {^ 2 ,m 2 }, {z/ 3 ,m 3 }}] = 
1 /■ /■ d d k l d d k 
ix d J J [k\ - m\\^ [kl - ml]"* [k\ - m 2 ]^ ' 

TJl[d,0, {{u u mi}, {z/ 2 ,m 2 }, {z/ 3 ,m 3 }}] = 

^) =1// (8) 

*™ 7r d 7 7 [A; 2 - m 2 ]-i [A; 2 - m 2 ]^ \k\ - m 2 ^ { 1 



and 



TBl[d, p 2 , {{^ 1?mi }, {^ 2 ,m 2 }}] = 

£(«*) = J_ / ^ f 9 ) 

tt<V 2 7 [A; 2 - m 2 ]-i - m 2 ]^ ' 1 J 

TAI^O,-^,!^}}] = 

Input for TARCER should be prepared in terms of the TFI-notation. 



4 



4 Elimination of numerators 



The TARCER-function TarcerRecurse is the principal function which, when 
applied to an expression involving the above integrals, performs the complete 
reduction to the set of basic integrals. In the first step the numerator of the in- 
tegrand in the TFI-integrals is simplified as far as possible by standard manip- 
ulations until an irreducible numerator of the form (Aki) a (Ak 2 ) b (pki) r {pk2) s 
results. The remaining integrals are of the form 

I [ 2s = TFl[d, p 2 , Ap, {a, b}, {0, 0, r, s, 0}, {{v u m,}, {u 5 , m 5 }}] . (11) 



The critical observation of Tarasov is that integrals containing irreducible 
numerators may be rewritten in terms of scalar integrals in higher space-time 
dimensions and that those can later be reduced again to scalar integrals in the 
original space-time dimension. Integrals of the form (11) we rewrite as 

I% rs = T^ s (p 2 ,Ap,{d},d + ) I$ 00 , (12) 



by employing an operator T that is a polynomial in the operator d + represent- 
ing a shift d — > d + 2 in dimension and in the mass derivatives dj = d/dm 2 . 
With the lightlike vector A present, the required T-operator is given by a 
generalization of eq. (25) of [lJQ 



T abrs (p 2 ,Ap,{d},d 
x exp 



d Yf d \ b ( d V 



d 



\id~ii) \id~i2) \id@i J \id(3 2 / 
i {(/V + 7 i(Ap)) Qi + (fop 2 + 72 (Ap)) Q 2 
+A(AP 2 + 2 7l (Ap)) Q n + (3 2 {fop 2 + 2 72 (Ap)) Q 22 
+ ((3i(3 2 p 2 + (/3i 72 + /3 2 7i)(Ap)) Qn} p 



(13) 



where the Q's are polynomials in dj 



Ai Qu = d 2 + d 4 + d 5 , -Q x = d 3 d 5 + d 4 d 5 + d 2 d 3 + d 3 d 4 , 
Ai Q 22 = di + d 3 + d 5 , -Q 2 = <9 4 <9 5 + <9 3 <9 5 + <9i<9 4 + d 3 d 4 , 

2iQ 12 = d 5 . (14) 



4 We have p = — d + instead of — 1/tt 2 d + due to the factor l/ir d in the normalization 
of the integrals in (12). 
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5 Recurrence relations 



Once all irreducible numerators are eliminated the next step taken by the 
function TarcerRecurse is to repeatedly apply the recurrence relations that 
reduce the exponents of the scalar propagators in the integrals until no further 
reduction is possible. All recurrence relations explicitly or implicitly given by 
Tarasov are implemented. 

For certain classes of on-shell integrals we found it necessary to supplement the 
recurrence relations by some additions in order to achieve maximal reduction. 
Those additional relations are listed in the following. The relations obtained 
by permutations of indices and masses are also implemented. The operators 
l 1 * 1 , 2 ± 1 . . . act on an integral by in-/decreasing the first, second, etc. index by 
one unit. 

For mi = and v 2 > 1: 

2 m\ (1 - v 2 ) v 2 2+ = 

(d - 2 v 2 ) (1 - v 2 ) 4% V3 - v x (-d + 2 v x + 2) l + 2~ 4%„ 3 . (15) 

For p 2 = ml, m 2 = m 3 = and v\ > 0: 



v x (d-2u 2 -2u 3 ) (d-v 2 -v 3 -l) 1 + 41L 3 = 

- v 2 (-d + 2 v 2 + 2) (-2 d + v x + 2 v 2 + 2 v z + 2)2 + 4% U3 . (16) 

For mi = m 2 = 0, p 2 = m\ and v\ > 1: 



(d - v_i - v 2 ) (-2 - d + 2 v x + 2 v 2 ) 
2m 3 2 + 
(3 d - 2 v x - 2 v 2 - 2 u 3 ) (-l-d + u 1 + u 2 + u 3 ) 
(2 d - 2 v x - 2 v 2 - v 3 ) (1 + 2 d - 2 v 1 - 2 z/ 2 - v 3 ) 

For mi = 0, m 2 = m 3 , p 2 = and z/i > 0: 



j(d) 

V\V2V3 



X 



1-41,3 • (17) 



- tt I/li/ 2 V3 



(d - 2 (-1 + i/x + f 2 )) (1 + d - v x - v 2 - v 3 ) 



2m 2 2 [d-2v x ) (l + d-2z/x 

(d-2 (-l + Z/i + Z/g)) 



x 



(2 + d - 2 z/i - z/ 2 - z/ 3 ) 
For m 1 = m 2 = m 3 , p 2 = m\ and V\ > 0: 



v 2 - v 3 ) 
K (d) 
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16 m\ v x {d-v x -v 2 - v 3 ) l + 4i 2Ki = 

_(l + 3 rf _ 3z/l _4z/ 2 )z/3l-3 + jS 2 , 3 
- v 2 (1 + 3 d - 3 v x - 4 v 3 ) 



+ (_i + 2 d - v x - 2 v 2 ) v 3 2S + 

- (-6 d 2 + 22 d i/i - 16 z/ 2 - i/ 2 + 6 d v 2 - 13 v x v 2 

- v 3 + 6 d v 3 - 13 vi v 3 - 4 v 2 v 3 ) J^l 2U3 

+ u 2 (-l + 2d-u 1 -2 u 3 ) 2+3- M . (ig) 



Note that in the last case all J^f U2U3 are eventually reduced to and (A^) 2 . 
That is, the number of basic two-loop integrals is smaller than in the general 
case. 

Other special configurations may require still further additional relations. 



6 Verification of recurrence relations 



As stated in [1] the derivation 'of these relations is rather tedious and for 
brevity of the presentation will be omitted'. Therefore, we rederived most of 
the recurrence relations. Furthermore, an automatic verification of recurrence 
relations for the integrals of type J or V with three and four propagators in 
the general mass case is possible via the kernel of the respective Mellin-Barnes 
type integral-representations. 

By replacing the massive scalar propagators in the momentum space integral 
with their respective Mellin-Barnes representation and interchanging integra- 
tions, as already described in [9] for the case ^ = 1, we obtain the Mellin- 
Barnes representation of J for general v i > 0: 

+ioo 3 

4?u = ?2-p /// «U(*i> *2, *o (- P 2 ) d — s n k 2 )" ds * > (2°) 

-too 1-1 



with the kernel 



^U( s i> s 2,s 3 )= (21) 

T(u + s-d) * T(- Si )T(d/2 -Uj- Si ) 
{ ' T(3d/2 - v - s) l\ r(^) 

where v = v x + v 2 + v 3 and s = si + s 2 + S3. As usual the integration contours 
separate the poles of the gamma functions in the left and right half planes. 
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The Mellin-Barnes representation of V contains four s-integrations in an anal- 
ogous manner with the kernel 

^ifu 2 u,u i (si,S 2 ,S 3 ,S 4 ) = (22) 
( l) 1 + u ^( l>1 + V -l + Sl + 33 ~ ^/ 2 ) r (^ - ^1 - ^3 - ^4 - gl - ^3 - S 4 ) 

+ ^3 + ^4 + si + s 3 + s 4 - d/2)T(d — i/i — i/ 3 — si — s 3 ) 

i> + s - d) r(-s 4 )i> 4 + s 4 ) ' r(- Si )r(d/2 - v A - sj 
x r(3rf/2 - 1/ - s) r(i/ 4 ) r(^) 

where z/ = v\ + z/ 2 + v 3 + z/ 4 and s = si + s 2 + S3 + s 4 . 

Insertion of the above integral representation for J into a vanishing combina- 
tion of J's as provided by a conjectured recurrence relation yields under the 
integral a sum of the form 

E 0n 1 n 2 n 3 (Sl,S 2 ,S 3 )(-p 2 )^ S K) Sl+ni (^2) S2+n2 (^3) S3+n3 ) (23) 
ni,n 2 ,n 3 

where the rij's run over a finite set of integers and n = ni + n 2 + n 3 . Here each 
represents a sum over ^'s with various values of the indices v { and d. By an 
appropriate shift of integration contours for each term separately accompanied 
by a change of variables Sj — > s« — one extracts a global factor: 

(V)^K) S1 (^2) S2 (^3) S3 E 0n 1 n 2 n 3 (si-ni,S2-n 2 ,S3-n3).(24) 

ni,n 2 ,n 3 



The sum in (24) uniquely determines the coefficients of the large p 2 expan- 
sion of the Mellin-Barnes integral over (24) in terms of (fractional) powers 
of m 2 /(—p 2 ). Therefore, if the recurrence relation is indeed correct and this 
integral vanishes then the above sum has to vanish identically. 

This observation provides a method to verify recurrence relations for both J 
and V, where a similar Ansatz applies. The method is suitable for recurrence 
relations that are valid for general values of m 2 and p 2 . It can also be adapted 
to the case when some of the m 2 or p 2 are zero. It is, however, not applicable 
to relations that hold only for particular parameter configurations, e.g., when 
some kinematical determinants vanish. 

This verification-procedure is implemented in TARCER through the functions 
CheckTJIRecurs ion and CheckTVIRecursion. The notebook TARCER. nb also 
contains comments for each particular recursion as to whether this check is 
applicable or not. 

Despite the complexity of the equations presented in [1] we found only one 



8 



non-obviousQ misprint: In eq. (67) the global sign of the right hand side has 
to be reversed. 



7 Usage and examples 



The Mathematical .0 notebook TARCER.nb contains the complete source code. 
Upon evaluation the recurrence relations are incorporated into the function 
TarcerRecurse. All T-operators (13) up to {a + b, r + s} < $RankLimit are 
constructed and the corresponding relations (12) explicitly generated. The 
value of $RankLimit can be set in the prologue section. 

The evaluation of TARCER.nb generates a binary file tarcer .mx containing all 
functions and definitions in internal Mathematica format. An explicit integral 
reduction only requires the latter to be loaded directly into Mathematica 3.0. 
A set of preproduced tarcer*. mx files for different operating systems is also 
available. Therefore one only needs to run the TARCER.nb notebook if one 
wants to increase $RankLimit or wishes to modify the program. 

TarcerRecurse is the main function that performs the integral reduction. Its 
usage is simply as follows: It can be applied to any expression involving the 
functions TFI, TVI, TJI, TBI and TAI whose first argument is a symbol d. 

For a simple example, consider the following on-shell integral where p 2 = M 2 : 

1_ f f d%d d k 2 

ir d J J k\ \k\ - M 2 ] [(k x - p) 2 - M 2 ] (k 2 - p) 2 [(h - k 2 ) 2 — M 2 ] ' 



Application of TarcerRecurse to the corresponding input form (4) yields: 



TarcerRecurse[TFl[d,M 2 ,{l, {1,M}, {1,M}, 1, {1,M}}]] 

3 (d-2) 2 (5d-18) (Ag M} ) 2 
32 (d - 4) 2 (d - 3) M 6 + 

(3d - 10) (3d -8) J{f [Af } i{ i [ o} i{ i [0 } 
8 (d - 4) (2 d - 7) M 4 

(3d — 10) (3d — 8) J{i^M},{l,AQ,{l,M} 

16 (d-4) 2 M 4 1 j 

5 The lower limit of the sum in (46) should read j = 1 and in (66) needs to be 
replaced by m 2 . 
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Note that the display format of TARCER explicitly shows the masses together 
with the indices. 

The function TarcerExpand inserts explicit results for some basis integrals as 
specified by the option TarcerReduce. A second argument to TarcerExpand 
must be given in form of a rule, like d — > 4 + e. Then an expansion of the first 
argument of TarcerExpand in the sole variable specified, here e, around will 
be performed. 

Applying TarcerExpand to % i.e. to the previous expression, yields: 



TarcerExpand[%, d — > 4 + e] 

{M 2 y- l S 2 E { 6 C(2) log(2) - (27) 



Here S e = e 7B ^ 4 <*' 2 is used.Q This result is of course well known, see [10] 
where the expansion up to 0(e 3 ) is given. 



As a somewhat more demanding example consider another on-shell integral 
with p 2 = M 2 containing Ak% in the numerator: 



1 r f d A k x d d k 2 (Ah 



dt,_ ( \U.\m 

(2f 



n d J J [kj - M 2 } 2 k\ (h - p) 2 [(k 2 - p) 2 - M 2 ] [(h -k 2 ) 2 - M 2 } 



For general m one can, with some effort, express this integral in terms of a 
moment integral, Laurent-expanded around d = 4: 



(MY" 6 (Ap) m S 2 J dx + | C(2) - 2(rf 1 _ 4)2 ) 6(1 -x) + 

log(l - x) _ 1 , _ \og(x) log(x) 

Uj+ 2(l-x) 4 l0gli X) 2(l + x) 2(l-x) 2 + 



(l - (i_ x f ) (2 l0g ^ l0g ^ ~~ ^ + 4 log2 ^)) ~ 

l - C(2) + log(x) log(l + x) + Li 2 (-oO) + 



1 - 

(1- 

n2 ) Li 2 (l-ar)) +Q(d-4)) . (29) 



l-x)VV2 
1 

2 (1 - a;} 



7s is the Euler-constant 
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TARCER was originally written to check this kind of general expressions for 
individual moments. Consider for example the moment m = 1 and apply 
TarcerRecurse to the input form (3): 



TarcerRecurse[TFl[d,M 2 ,Ap, {1, 0}, {{2,M}, 1, 1, {1, M}, {1, M}}]] 

(d - 2) 2 (27 d 4 - 404 d 3 + 2175 d 2 - 4902 d + 3776) Ap (A ( ^ M} )' 



(d - 2) (3 d - 10) (3 d - 8) Ap jf? 



256 (d-5) 2 (d-4) 2 (d-3) M 8 

{1,M},{1,0},{1,0} 



32 (d - 5) (d - 4) (2d - 7) M 6 
d (3 d - 10) (3 d - 8) Ap jff M}i{1 , M}i{1 , M} 



+ 



(30) 



128 (d - 4) 2 M 6 
Expansion with TarcerExpand yields: 

TarcerExpand[°/„, d — > 4 + e] 

(M 2 )- 2 Ap S 2 ( 3 log(2) C (2) - ^ --^-^3)+^ (31) 

Which agrees with the analytic integration of the moment-integral above for 
m — 1. 



For a massive off-shell example consider the integral 

Iff d%d% (kj) 2 

W J [k 2 - m\\ [(h - p) 2 - m 2 } [(h - k 2 ) 2 - m 2 ] ' ^ ] 



denoted Y^ A in [8]. Acting with TarcerRecurse on the input form (2) yields: 
TarcerRecurse [TFl[d, p 2 , {2, 0, 0, 0, 0}, {0, {l,m 4 }, {l,m 2 }, 0, {l, m 3 }}]] 

(d) a (d) 

L {l,m 2 } A {l,m 3 } 



(d (p 2 + m 2 + ml + 7 mi) - 12 m 2 ) A«> A& 



3 (3 d - 4) 



+ 



(d (p 2 + m\ + 7 m\ + mj) - 12 mj) Ag ma} 
3 (3 d - 4) 

((7 d - 12) m 2 + d {p 2 + m 2 + m 2 )) Ag ms} A^ 
3 (3 d - 4) 
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-(((d - 2) d p 4 + 2 (d - 4) (d - 3) (m 2 + m 2 ) p 2 + 



3 (d-2) (3 d-4) 

(d (25 d - 102) + 96) m\ - 2 (d - 3) (d m 4 + 

2 (7 d - 12) m 2 m 2 + d m 4 ) + 2 m 2 ((19 (d - 4) d + 72) p 2 - 

(d - 3) (5 d - 12) {ml + m 2 ))) jg^ HW{W ) + 

' '2 m 2 (-3 (3 d-4) m 4 + 2 ((18 - 11 d) p 2 + 



3 (d- 2) (3 d-4) 

2 (d - 3) ml) m 2 2 - 12 m\ (-p 2 + 3 m 2 + m 2 ) + 

d (_ p 4 _ 4 m 2 p 2 + m 4 + 9 m 4 + 22 m 2 m 2 )) j;? m3}{lim4}{lim2} ) 

' 2 ml (-3 (3 d-4) m 4 + 2 ((18 - 11 d) p 2 + 



3 (d- 2) (3 d-4) 

2 (d - 3) m 2 ) m 2 2 + 3 (3 d - 4) m\ + 2 m 2 ((11 d - 18) m 2 - 

2 (d - 3) P 2 ) + d (m 4 - p 4 )) jg^HLma}!!^}) - 

' (8m 2 [m 2 2 -p 2 ) ((2d-3)p 2 + (2d-3) m 2 



3 (d-2) (3 d-4) 

{2,m2}{l,rri4}{l,m3} 



(d - 3) (m 2 + ml)) jg m2H1 , m4H1 , m3> ) (33) 



which is a much shorter expression than the one given by eq. (B.2) in the 
appendix of [8]. Moreover, no spurious massless propagator is introduced. 

These and further examples can also be found in the material on the WWW- 
sites listed in the program summary. 



References 



[1] O.V. Tarasov, Nucl. Phys B502 (1997) 455-482; jhep-ph/97033l9| . 

[2] O.V. Tarasov, Phys. Rev. D54 (1996) 6479; |hep-th/96060Tl . 

[3] J. Lascoux and T. Regge, In: Battelle Rencontres 1967, eds. CM. DeWitt and 
J. A. Wheeler (W.A. Benjamin, New York, 1968). 

[4] R. Hamberg and W.L. van Neerven, Nucl. Phys. B379 (1992) 143. 
R. Hamberg, University of Leiden, thesis (1991). 

[5] S.A. Larin, F.V. Tkachev and J. A.M. Vermaseren, NIKHEF-H-91-18 (1991). 
S.G. Gorishnii, S.A. Larin, L.R. Surguladze and F.V. Tkachev, Comput. Phys. 
Commun. 55 (1989) 381-408. 

[6] J. Fleischer, O.V. Tarasov, Comput. Phys. Commun. 71 (1992) 193-205. 

[7] G. Weiglein, R. Mertig, R. Scharf and M. Bohm, In: New computing techniques 
in physics research II, ed. D.Perret-Gallix (World Scientific, Singapore, 1992) 
p.617. 



12 



[8] G. Weiglein, R. Scharf and M. B6hm, Nucl.Phys. B416 (1994) 606 
ph/931035£ 

[9] F. Berends, M. Buza, M. Bohm, R. Scharf, Z. Phys. C 63 (1994) 227. 
[10] D.J. Broadhurst, Z. Phys. C 54 (1992) 599. 



13 



